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In this work we study the presence of kinks in models described by a single real scalar field in 
bidimensional spacetime. We work within the first-order framework, and we show how to write first- 
order differential equations that solve the equations of motion. The first-order equations strongly 
simplify the study of linear stability, which is implemented on general grounds. They also lead to a 
direct investigation of twinlike theories, which is used to introduce a family of models that support 
the same defect structure, with the very same energy density and linear stability. 
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I. INTRODUCTION 

Kinks and lumps are defect structures of current in- 
terest in high energy physics [l], 01 and in other areas of 
nonlinear science l3l|. In high energy physics, in partic- 
ular, kinks and lumps appear in models described by a 
single real scalar field in bidimensional spacetime, and 
they are static solutions of the corresponding equations 
of motion. Recent examples of models where kinks play 
important role can be found, for instance, in Refs. 

See also Ref. [1^, where one studies bell-shaped or 
lump-like solutions in models described by a single real 
scalar field. 

In Cosmology, defect structures are usually related to 
the existence of phase transitions in the early Universe. 
Motivated by Cosmology, however, another direction of 
study concerns models where kinematic modifications of 
the scalar field are introduced aimed to contribute 
to explain the present accelerated expansion of the Uni- 
verse. Recent investigations of scalar fields with modified 
kinematics have been presented in [l7l - [23 | with distinct 
motivations, in particular to search for the presence of 
topological structures which live in a compact region, 
and for applications to the braneworld context, with the 
five dimensional spacetime engendering a single extra di- 
mension of infinite extent. For other recent applications 
within the brane-world context, see, e.g., Ref. |2a |. 

In the study of models with modified kinematics, in 
[26l | the authors noted that it is possible for k-defects 
to masquerade as canonical scalar field solutions. That 
is, given a standard scalar field model, the topological 
defect profile and corresponding energy density can also 
appear in a k-field theory. In this case the two models 
are twins of each other, and the idea was further explored 
very recently in [13 , and also in [28|. Motivated by these 
recent results on twinlike models '26'-'28'], in the present 
paper we introduce and study a new family of models, 
described by a single real scalar field that can accommo- 
date interesting results. Before doing that, however, in 
the next Sec. |ll] we briefly review some results obtained 
in the case of standard models. We then move on to the 
main results of the present work, which are described in 



Sec. Illli where we discuss the new family of models, and 
in Sec. IIVI where we study twinlike models. We end the 
paper in Sec. |Vl with some comments and conclusions. 



II. STANDARD MODELS 

We start with a scalar field cj), described by the stan- 
dard Lagrange density 



(1) 



where U{(j>) is the potential, which is used to identify the 
particular model under consideration. In this work we 
deal with bidimensional spacetime, with metric such that 
x'^ = xo = t and x^ = —xi = x. For simplicity, we use 
dimensionless field and coordinates, so all the quantities 
are also dimensionless. 

We suppose that the above model engenders static so- 
lution, with (j) = 4>{x) that obeys the equation of motion 



(2) 



where prime stands for the derivative with respect to x, 
that is, 4>' = d4>/dx, and we are using t/^ to represent 
dU/dcf). The energy density of the static solution is given 

by 



six) = + U{<t>) 



(3) 



We can integrate the equation of motion ([2]) to obtain 

= 2C/(0) + C (4) 

where C is a constant. If one considers static field with 
finite energy, we have to take C = 0. This leads to the 
equation 



^'2 = 2[/(0) 



(5) 



which is a first-order differential equation that describes 
static solution with vanishing stress component of the 
energy-momentum tensor. We can use ([S]) to introduce 
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an interesting result, concerning the equipartition be- 
tween the gradient and potential portions of the energy 
density. We take 

= \<i^'^ (6) 

as the gradient part of the energy density, and 

= C/(0(x)) (7) 

as the potential part. We see that 

e{x) = £g{x) + £p{x) = 2eg{x) = 2£p{x) (8) 

for the static solution that solves the first-order differen- 
tial Eq. ([S|). Similar expressions also works for the energy, 
and we can write E = Eg + Ep = 2Eg = 2Ep. 

This result is general, and it is valid under the above 
kinematics, which represents the standard situation, de- 
scribed by the Lagrange density ((TJ. However, if one 
changes kinematics, the energy equipartition is in gen- 
eral lost. In this work we focus our attention on a new 
family of models, which represents an interesting possi- 
bility to depart from the standard case, since it will lead 
to new results on twinlike field theories. 



III. NEW MODELS 

Before introducing the family of models, let us intro- 
duce the quantity 



(9) 



Thus, we can write the standard Lagrange density ([T]) in 
the form 



C = X- U{<j>) 



(10) 



Usually, one goes beyond the standard case with the 
two possibilities: we take F{X) as a general function of 
X and we consider the case 



C = F{X) - U{^) 



(11) 



Here we note that to specify the model, we need to have 
both U{(f) and the function F{Y), and we see that in the 
case 



F{Y) = 1 + Y 



(15) 



we get back to the standard model. 

The equation of motion for the family of models has 
the form 



{FYd^4>) + {F- YFy) C/0 = 

where Fy = dF/dY. 

The energy-momentum tensor is given by 



-g^,U{^)F{Y) 



or, in components 

Too 

Toi 
Til 



Fy4)^ + U{(f>)F 

FY(t>'4> 

Fy(I)'^ - U{(j))F 



(16) 

(17) 

(18a) 
(18b) 
(18c) 



where we are using dot to represent time derivative. Let 
us now suppose that {vi, « = 1, 2, n} is a set of static 
and uniform solutions of the equation of motion. This 
means that U^{vi) has to vanish. Also, from the energy 
density we take U{vi) — to make the energy of the static 
and uniform solutions vanish, in the same way they do 
in the standard case. 

Since we are dealing with a new family of models, we 
guide ourselves with the null energy condition (NEC), 
that is, we impose that T^^n^n'^ > 0, where is a null 
vector, obeying gf_i,^n'^n^ — 0. This condition leads to 
following restriction Fy > 0, for the general field config- 
uration (f>{x, t) which solves the equation of motion (jl6l) . 

Let us now search for defects structures, considering 
the case of a static field configuration, cf) = (f>{x). In this 
case, the equation of motion ([T5|) changes to 



{2FyyY + Fy) 



uAf 



FyY + 2FyyY^) 



This equation can be integrated to give 



2YFy - F ■ 



C 



C/(0) 



(19) 



(20) 



where C is a constant. In the case of static solution, Y 
becomes 



C = U{^)F{X) 



(12) 



which includes the Born-Infeld generalization. In this pa- 
per, however, we choose a different possibility, described 

by 



where 



C = -U{^)F{Y) 



U{(b) 2 U{(b) 



(13) 



(14) 



Y 



2{/(0) 



(21) 



Thus, the above Eq. ([^0]) is a first-order equation, and we 
see that its solutions make Tn = C. Thus, the constant 
of integration that leads us to the first-order equation 
(1^01) also identifies the stress component of the energy- 
momentum tensor. 

The first-order equation ((20)) can be written as follows 



C/(0) 



(22) 
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where the function G is such that G-'^{Y) = 2YFy - F. 
Note that when C vanishes, we get to the case of stressless 
solutions, and we obtain an algebraic equation for Y, such 
that F = 2YFy; for G(0) — c, with c constant, real, this 
leads to a simpler relation 



= 2cf/((/)) 



(23) 



and so we have Y = c. Note that the solution of this 
equation is the solution 0^ (x) of the equation of motion 
(O of the standard model, with the position changed as 
X — > -y/c X, that is, we can write 



<p{x) = (psiVcx) 

The thickness of the solution is given by 
6 - SJV^ 



(24) 



(25) 



so the solution is thicker or thinner depending on the 
value of c being lesser or greater than unit. We note that 
the case of c negative is not allowed; also, only stressless 
solutions have the form given by Eq. 

We see that the profile of the stressless solution is sim- 
ilar to the profile of the solution in the standard model. 
Also, the energy density of the static solutions have the 
form Too = U{(j))F{Y). In the case of stressless solutions 
we have C = and we get 



Too = F{c)U{(t>) 

Here we can write 

E^F(c) / U{(l){x))dx 



or better 



E 



Fjc) 



U[Uv))dy^^E, 



(26) 



(27) 



(28) 



where Eg is the energy of the standard solution. 

The solutions with a non- vanishing Tn are completely 
distinct from the corresponding solutions of the standard 
model, since they do not have the form given by Eq. ([M]) . 
Here we recall that for Tn — C, only the stressless so- 
lutions are stable [l^. Usually, the energy of the other 
solutions are divergent, and the solutions have oscillatory 
or divergent profiles as they do, for instance, in the stan- 
dard model. However, it is interesting to observe that 
for some choices of the function F{Y) we get to the case 
of Til = C ^ 0, with solutions with topological behavior 
and finite energy. We further illustrate this new possibil- 
ity with an example below. 



A. Linear Stabilty 

Let us now investigate linear stability of the static so- 
lution. We introduce small fluctuations r]{x,t) about the 
solution in the usual way: we write (j){x, t) = 4>{x) + 



r]{x,t), where (/)(a;) represents the static solution. The 
equation of motion (jl6p allows obtaining, up to the first- 
order power in 77, 



d.„ FYd^v-^d^ci)d^<i>d''v 

U{(p) 

" {FyY-F)yY 



U^^ (FyY - F) 



rj 



(29) 



Since (j) — (j){x) is static solution, we take •q{t, x) 
'q{x) cos(cj t) to obtain 

-[{2FyyY + Fy)v'{x)-^ = 

ul 



{FyY-F) 



( U^FyyY^' V 
V UU) J 



{FyY - F)yY 



lo^Fy 



r]ix) 



In the case of the standard model we have F{Y) 
and so we get 



(30) 



- Vs{x) + C/00|0=^^(^)??s(a;) = (31) 

which describes the standard expected. 

In the general situation, for stressless solutions we use 
the equation (^0)1 with C = to obtain 



T]"{x) + c t^0^|^=^^(^^)?7(a;) - Vix) 



where 



A' 



2FyyY + Fy 



(32) 



(33) 



To ensure hyperbolicity of the differential equation, we 
have to impose the constraint that is constant and 
positive; see, e.g., Ref. [l^, and remember that A^^ is 
the sound speed. 

If we introduce the transformation z — ^/ca;, we can 
rewrite the above equation (j32p in the form 



'I0=0.(^) 



r^{z) (34) 



which has the very same form of the Eq. (PT|) , which de- 
scribes stability in the standard model. Thus, we can find 
a one-to-one correspondence between the eigenvalues (w) 
and eigenstates (77) of the general theory and the eigen- 
values {lOs) and eigenstates {rjs) of the standard model. 
They are related by 



LO = UJgA^/c 

ri{x) = ris{\fcx) 



(35a) 
(35b) 



This is an important result, since the behavior of the 
stressless solution in the general model given by the La- 
grange density (fT3)) is similar to the standard model, dif- 
fering only by a scaling in x and in the eigenvalues. In the 
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standard model, if we choose a potential U{(j)) that leads 
to stable static solution, then the general theory also sup- 
ports stable static solution, under the assumptions that 
^2 > and c> 0. 

To end the study of stability, let us go further and in- 
vestigate the second-order contribution in rj which appear 
in the energy density. They are given by 



T, 



1 



LoTj Fy sin [ut) 



+ -FyA2(,7'2 + cC/00r/2)cos2(wt) (36) 

Since the contribution is time-dependent, we take the 
average in time to get 



Fy 



Now, using the Eq. p2p we obtain 



2 2 

uj rj 



(37) 



(38) 



We integrate this expression to get 



<£:(2)> = FyU)'^ / T]'^{x)dx 

J — oo 

/oo 
rjl{z)dz 
-OO 



F{c) 



Using the relation (fT6|) , we can write 



<eP> 



(39) 



(40) 



which shows that for an admissible F{Y), it is easier or 
harder to perturb the solution, depending on A being 
lesser or greater then unit, respectively. 



B. Examples 

Let us now specify the function F{Y) in order to illus- 
trate how the formalism introduced above works. First 
we recall that for F{Y) = 1 + F, we get back to the 
standard model. 

We introduce a model that can be understood as an ex- 
tension of the standard model where the dynamics engen- 
ders a nontrivial factor, of the form £ = g^cf), X)X—U{(j)). 
In particular, we choose the following Lagrange density 



aU{cj)) (41) 



Here we suppose that rt > 1 is real parameter, and a and 
/3 are real and positive. Note that for n = 1, the model 
describes standard kinematics. This model is different 







FIG. 1: Plot of the solutions in the {<f>, <jf>') plane, for C = 
(black, solid line), C < (blue, dashed-dotted line) and C > 
(red, dashed line), from left to right and up to down, for 
n = 1, 1.1, 1.5, and 2, and for a = [3 = 1. 



from the model introduced in Ref. [T^ and also inves- 
tigated in Ref. [U [l^. It can be written in the form 
given by Eq. if we identify F{Y) as 



F{Y)=a + l3Y\YY 



(42) 



In this case, both the NEC and hyperbolicity conditions 
are satisfied. Using the first-order Eq. (I^H)) . we get 



•=2 



1 



f3{2n~ 1) 



aU{(t>Y 



(43) 



We note that for n> 1, we can find topological solutions 
even though we consider C > 0. Such solutions have the 
same kink profile of the stressless solution, but now they 
are thinner than the stressless solution. This behavior 
reflects the fact that solutions with C greater then zero 
must have energy greater then the stressless solution, so 
they must be thinner than the stressless solution. In 
Fig. 1 we depict the above expression P5)) in the (0, </>') 
plane for the standard potential 



(44) 



We take C = in the above Eq. (|43D to consider the 
stressless solution. We get 



(45) 



/3(2n- 1) 
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and so 



E 

'K 



2n-l 



n 



(46) 



We can fix the thickness of the defect as in the standard 
model (dl), choosing c — 1; this leads us to the choice 
a = {2n — l)/3. The energy behaves as 



E 

e: 



2n-l 



(47) 



and so it decreases with increasing n. The stability of 
the solution does not depend on a and /3, and we get 



(48) 



Using the relation (PU)) . we see that when n increases, it 
becomes harder and harder to destabilize the solution. 
We introduce another model, described by 



(49) 



If we make the transformation = 
[/(0(T))a^T5^T, we get C = -U{q^iT))^l - d^Td^^T 
which can be used to describe the dynamics of tackyonic 
brane; see, e.g., Ref. (29l |. 

The Lagrange density (H^ has the form of the general 
model ([13)) . where the F{Y) function is 



F{Y) = vT+W 



The condition ([2^ gives 



Q2 



C/(0) 



(50) 



(51) 



The model has no solutions when C is nonzero. For a 
stressless configuration, C = 0, we see that Y ^ oo. The 
solution has the step-like behavior 



(x) = Wl + (V2 - vi)H{x) 



(52) 



where vi and V2 are minima states, and H{x) is the Heav- 
iside function. It has the profile of a topological structure 
with vanishing thickness. The energy density correspond- 
ing to (l52|) has the form Too = iV2/2)S{x — xq), where 
xq is the center of the solution, and the energy is finite, 
given hy E = V^/2. 

In this case we can calculate that 



A' 



1 



1 + 2Y 



(53) 



which vanishes for y — > oo, that is, for the stressless 
solution, with C — 0. This is similar to what happens 
in the cuscuton model [M], the difference is that here 
this happens with the stressless solution, while in the 
cuscuton model every solution presents this behavior. 



In order to regularize the solution we follow the case 
considered in 30] and we add a constant to F{Y), that 
is, we rewrite (|50p as 



F{Y) = b VT+2Y- - 



(54) 



where a and b are real parameters. The solution now 
obeys 



Y = 



(55) 



and the model supports valid solutions for a > 1, with 
energy 



E = ^y/a^ -1 
V2a 



(56) 



We note that obeys the same expression given by Eq. 
((53| : therefore, — a^^, and a is the sound speed. 



IV. TWINLIKE MODELS 



According to some recent results j26l - l28l |. two distinct 
models Ci and £2 are twinlike models if they support 
static solutions with the same profile, and with the very 
same energy density. 

The family of models here introduced is controlled by 
F{Y) and U{<j)). Thus, let us now introduce conditions on 
the function F{Y) for a theory given by the action 
to be twin of the standard theory, where F{Y) — 1 + Y. 
Comparing the former Eqs. ([5]) and (j23l) . for c = 1 we see 
that they have the same static solution, and this fulfills 
one of the two condictions for the models to be twins. In 
this case, the value of Y at the static solution has to equal 
unit, that is, we have to impose that y = 1 for the static 
solution. Next, we compare Eqs. dH), (O, and (^5)) to get 
to F{1) — 2, in order for the static solutions to have the 
same energy density. From this relation, we obtain that 
Fy(Y — 1) should also be unit, that is, ^V(l) = 1- 

These are the general restrictions on F{Y) to make the 
model twin of the standard model. For instance, in the 
first example studied before, if one writes 



F(y)^^(i 



y|y|' 



2n- 1 



(57) 



one gets to the Lagrange density 

n-l 



C 



2"n 



2n- 1 



[/(0) (58) 



which is twin of the standard model with the same U{(f)). 
Similarly, we can rewrite the model with F{Y) given by 
Ea.([54|) in the form 



C/(0) 



-73 



(59) 



6 



with a and b chosen to offer another family of twin models 
to the standard model with the same U{4>). 
We go on, and we consider the function 

F{Y) = l + Y+^il-Yf (60) 

with a real constant. This is the simplest polynomial 
function which makes the model twin of the standard 
model, with F{Y) = 1 + Y. For this choice of F{Y), 
we note that = 1, so the two models are twins, and 
they also have the very same stability features. This 
result is interesting since it shows that there are twinlike 
models that cannot be distinguished by linear stability, 
differently from the systems investigated in Refs. 26, 27]. 

V. COMMENTS AND CONCLUSIONS 

In this work we studied the presence of kinks in models 
described by a single real scalar field in bidimensional 
spacetime. We introduced a family of models, which is 
defined by the two functions F{Y) and U{(p), and has 



the general structure C = —U{4>)F{Y)^ where we used 
Y = -a^0a''0/2f/((/)), with U{(j)) being a function of the 
real scalar field 4>. We note that for F(Y) = 1 + 1" the 
model becomes the standard model, described by Eq. ([T}. 

As an interesting result, we introduced the family of 
models defined by 

F{Y) = 1 + Y + ^{1-Yf 

which provides twinlike models to the standard model, 
with F{Y) = 1 + y. For such models, the presence of 
defects and their stability depend on the choice of U {(/)), 
but for each U{4>), the standard and twin models have 
the very same stability features. 

The issue concerning twinlike models deserves further 
investigations, and we are now considering the possibil- 
ity to extend the idea to the case of more sophisticated 
defects, such as vortices and monopoles, with the addi- 
tion of Abelian and non-Abelian gauge fields. Another 
issue concerns the presence of fermions, leading us to 
supersymmetric extensions of the twinlike models here 
introduced. 
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